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Abstract — We consider the design and analysis of the 
efflciently-encodable rate-compatible (E 2 RC) irregular LDPC 
codes proposed in previous work. In this work we introduce 
semi-structured _E 2 i?C-like codes and protograph E 2 RC codes. 
EXIT chart based methods are developed for the design of 
semi-structured £' 2 i?C-like codes that allow us to determine 
near-optimal degree distributions for the systematic part of the 
code while taking into account the structure of the deterministic 
parity part, thus resolving one of the open issues in the original 
construction. We develop a fast EXIT function computation 
method that does not rely on Monte-Carlo simulations and can 
be used in other scenarios as well. Our approach allows us to 
jointly optimize code performance across the range of rates under 
puncturing. We then consider protograph E 2 RC codes (that have 
a protograph representation) and propose rules for designing 
a family of rate-compatible punctured protographs with low 
thresholds. For both the semi-structured and protograph E 2 RC 
families we obtain codes whose gap to capacity is at most 0.3 
dB across the range of rates when the maximum variable node 
degree is twenty. 

Index Terms — E 2 RC codes, EXIT chart, semi-structured 
LDPC codes, capacity approaching, joint optimization, rate- 
compatible, puncturing performance, protograph LDPC codes, 
density evolution 



I. Introduction 

Low-density parity-check (LDPC) codes [1] have found 
widespread acceptance in different areas due to their superior 
performance and low complexity decoding. In this paper, we 
investigate rate-compatible punctured LDPC codes that have 
the flexibility of operating at different code rates while hav- 
ing a single encoder-decoder pair. Rate-compatible punctured 
codes are defined by specifying a systematic mother code that 
operates at the lowest code rate. The parity bits of higher 
rate codes in a rate-compatible code family are subsets of 
the parity bits of lower rate codes. A number of papers have 
investigated issues around the design of good rate-compatible 
punctured LDPC codes. The work of [2] presents methods for 
finding optimal degree distributions for puncturing. In [3] [4] 
[5], algorithms for finding good puncturing patterns for a given 
mother code were proposed. There have also been attempts 
to design mother codes (along with puncturing patterns) with 
good performance under puncturing [6] [7] [8]. 

E 2 RC codes introduced in [6] are linear-time encodable 
and have good puncturing performance across a wide range 
of code rates. In this work we present systematic approaches 
for the design and analysis of £' 2 _RC-like codes. Let H = 
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[H\\H2\ denote the parity check matrix of a systematic LDPC 
code where Hi denotes the systematic part and the parity 
part. We address the design of two types of codes in our work 
as explained below. 

i) Semi-structured E 2 RC-like codes. In these codes the 
parity part H2 is deterministic. We use the lower trian- 
gular form introduced in [6] and introduce a protograph 
structure for the H2 part. An example is shown in Fig. Q] 
We assume a random edge interleaver between systematic 
variable nodes and check nodes, which divides the code 
into a structured part and an unstructured part, as shown 
in Fig.Q] We solve the problem of finding optimal degree 
distributions for the unstructured part in this case for 
optimizing the rate-compatible codes at any specified 
punctured code rate(s). 

ii) Structured E 2 RC-like codes. These codes are protograph 
codes as introduced in [9]. The distinguishing feature 
is that the parity part of the protograph has an E 2 RC 
structure. We demonstrate that very good rate-compatible 
punctured code families can be obtained using the design 
rules we propose for the protograph construction. The 
protograph structure is especially valuable in practical 
applications as it allows parallelized decoding and re- 
quires significantly less storage space for the description 
of the parity-check matrix than unstructured codes when 
circulant permutations are used. 

We obtain semi-structured E 2 RC codes that have a small gap 
to capacity across the range of puncturing rates. Furthermore, 
we present optimized quasi-cyclic protograph codes based on 
the E 2 RC structure and demonstrate that very good perfor- 
mance can be obtained with them. 

This paper is organized as follows. In Section [II] we 
briefly discuss the main contributions of our work. Section 
|IV] presents our new method for the design of semi-structured 
E 2 RC codes. We also discuss the method of predicting the 
puncturing performance of semi-structured E 2 RC codes and 
the joint optimization of our codes at any specified punctured 
code rates. We explain the construction of protograph E 2 RC 
codes in Section |VJ and Section [VI] outlines our conclusions. 

II. Main Contributions 

We first outline the issues left unresolved in the work of [6]. 

a) The original construction of E 2 RC codes proposed the 
special H2 (parity part) structure of the parity-check matrix H. 
However the design of appropriate degree sequences for the 
H\ (information part) based on the constrained H2 structure, 
was not discussed. In [6], the authors used degree sequences 
designed for standard irregular codes and constructed H to 
match these distributions as closely as possible. 
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Fig. 1. Tanner graph representation of E 2 RC codes. 



b) The construction technique did not provide any means of 
optimizing code performance at any particular puncturing rate 
or across all rates simultaneously. 

c) As pointed out by an anonymous reviewer, the original 
E 2 RC codes suffer from high error floors [10] at the mother 
code rate. As shown in [10], this is because the H2 structure 
causes the maximum check node degree to be large. 

d) The original E 2 RC codes work with completely random 
interleaves, that are hard to implement in practice. 

In this paper, we resolve each of the issues discussed above. 
We briefly overview the main contributions below, 
i) Systematic design techniques for E 2 RC-like codes. 
Note that the analysis of E 2 RC codes does not follow 
directly from the analysis of related codes such as systematic 
IRA codes [11] [12]. This is because the structured part of 
IRA codes is symmetric while that of E 2 RC codes is quite 
asymmetric. In [12], four methods were proposed for the 
design of IRA codes. The first two methods implicitly assumed 
one edge type in the accumulator part which was justified 
by the symmetry of the part. Together with a one-parameter 
approximation of the message distribution function, Gaussian 
or BEC approximation, these two methods yielded almost 
closed-form equations of density evolution. However, one- 
edge type assumption turns out not accurate enough for the 
structured part of E 2 RC codes because of its asymmetry. In 
the latter two methods in [12], Monte Carlo simulations were 
used for generating the EXIT function of the structured part 
of IRA codes. The Monte Carlo simulation based method is 
accurate for computing EXIT functions of both symmetric 
and asymmetric constituent code components by taking the 
structure of the code component into account. When we design 
semi-structured E 2 RC codes using EXIT chart, we take into 
account the complete structure of the deterministic part of 
E 2 RC codes to compute the EXIT function as presented in 
Section [TV] Instead of resorting to Monte Carlo simulations, 
we propose a fast and analytical method for computing EXIT 
functions by solving a set of equations. We use multiple 
edge types [13] for the structured part of E 2 RC codes, one 
edge type for each edge in the protograph representation. So 
instead of having only three equations (equations (19) (20) 
(21) in [12]) from the structured part of IRA codes, we have 
2\Er\ + \El\ equations from the structured part of E 2 RC 
codes for density evolution. As demonstrated by simulations 



and the threshold predictions, this introduces a systematic 
method towards the design of semi-structured E 2 RC codes 
with better performance than the original E 2 RC codes. 

ii) A fast technique for EXIT function computation of code 
components based on protographs. 

Note that usually EXIT functions are computed via Monte- 
Carlo simulation, which tends to be time-consuming. In this 
work we present a general technique for computing EXIT 
functions of code components with a protograph structure. 
This greatly speeds up the code design process. While we 
applied it to the design of our semi-structured E 2 RC -like 
codes, it can be applied for any protograph like components, 
e.g. we can apply it to find the EXIT function of the H2 part of 
the IRA code by working with its protograph representation. 

iii) Simultaneous optimization of code performance across 
multiple rates. 

By exploring the E 2 RC structure and its designed puncturing 
pattern, we propose the design of good rate-compatible punc- 
tured codes so that the gap to capacity across the entire range 
of rates can be controlled. To the best of our knowledge, the 
current literature does not address this point. 

iv) Alleviating the high error floor problem of the original 
E 2 RC codes. 

In our design of semi-structured E 2 RC codes, we impose a 
protograph structure on the H2 part, which corresponds to the 
H2 part of a very short original E 2 RC code. This ensures that 
the maximum check node degree remains low, thus preventing 
the high error floors that occur in the original E 2 RC codes at 
mother code rate. For a related approach see [10]. 

v) Design of high-performance codes based on protographs. 
Codes with completely random interleavers are too complex 
from the point of view of implementation in hardware. In 
this work, we design protograph E 2 RC codes where both 
the Hi and the H2 parts have a protograph structure. We 
propose design rules for generating a family of rate-compatible 
protographs with good threshold properties at all punctured 
rates. Finally, we demonstrate codes with performance better 
than the original E 2 RC codes, that are obtained by replacing 
the protograph edges by circulant permutations. 

III. Background and Related Work 

An LDPC code can be defined by a parity-check matrix or 
equivalently by a bipartite (or Tanner) graph representation. 
For the bipartite graph representation, we follow the conven- 
tion that a blank circle represents an unpunctured variable 
node participating in the transmission and a filled circle 
represents a punctured variable node not participating in the 
transmission. The asymptotic threshold of LDPC codes can 
be found by performing density evolution [14] [15] [16] [17] 
on the degree distribution pair. However, for LDPC codes 
with structured components such as IRA codes and protograph 
LDPC codes [9], the density evolution analysis needs to take 
the underlying structure into account. This can be handled by 
classifying edges into different types [13] and also by using 
EXIT charts [18]. Protograph LDPC codes start with a small 
mini-graph (called a protograph) and construct the LDPC 
codes by replacing each edge in the protograph by a random 
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permutation of a fixed size. They can be considered as a 
subclass of the multi-edge type LDPC codes [13]. Fast density 
evolution based on the reciprocal channel approximation [19] 
can be performed on protographs to determine their asymptotic 
threshold. 

A. Efficiently Encodable Rate-Compatible LDPC Codes 

We now briefly overview the E 2 RC codes introduced in 
[6]. Let H — \H1\H2\ denote the parity-check matrix of an 
E 2 RC code in systematic form. We say that a parity node 
in H2 is fc-step recoverable (or k-SR) if it can be recovered 
in exactly k iterations of iterative decoding assuming that 
all the parity bits are punctured and all the systematic bits 
are known (Fig. [2] shows an example). Intuitively, a large 
number of low-SR nodes tend to reduce the required number 
of decoding iterations in the high SNR regime and result in 
good puncturing performance. 

1-SRnode 2-SR node 3-SR node 




Fig. 2. The figure shows an example of a 1-SR, 2-SR and 3-SR node. 

In [6], the submatrix H2 consists of exclusively degree-2 
and degree- 1 nodes. Moreover, when the number of parity 
nodes is a power of two, half the nodes in H2 are 1-SR, one- 
fourth are 2-SR and so on. The special structure of H2 for 
E 2 RC codes allows linear-time encoding and results in good 
puncturing performance with a puncturing pattern, where 1- 
SR nodes should be punctured first, 2-SR nodes be punctured 
next and so on depending upon the rate requirement. 

The E 2 RC codes have good puncturing performance at 
relatively short block lengths. However, when the block length 
gets large, the structure of H2 may induce a large spread in 
the check node degree distribution that may cause a loss of 
performance. In recent work, Song et al. [10] showed that 
E 2 RC codes exhibit high error floors at their mother code rate 
and claimed that this stems from their dispersive right degree 
distribution and high maximum right degree. They presented 
a modified approach that fixes the high error floor problem. 
In Section IIV-BI we show that our approach also effectively 
eliminates the high error floors of E 2 RC codes at the mother 
code rate. In fact we obtain codes whose performance is 
slightly better than those in [10]. 

B. EXIT Chart Overview 

EXIT charts [20] were first proposed for understanding 
the convergence behavior of iteratively decoded parallel con- 
catenated codes, and were later generalized to the analysis 
of LDPC codes [18] [21] [22] [23]. The components of an 
EXIT chart are the EXIT functions of the constituent code 
components of the iterative decoder, which relates the a priori 
mutual information available to a code component, denoted I a 
and the extrinsic mutual information generated after decoding, 
denoted Ie- The advantage of EXIT charts is that the code 



design problem can be reduced to a curve fitting problem 
between the code components (usually two in number). 

For log-domain belief propagation decoding of unstructured 
LDPC codes, if the incoming messages to a variable node v 
of degree d v are assumed to be Gaussian and independent, the 
EXIT function for the code component involving all variable 
nodes is given by [21] 

I E, V {I A, V, ^nch) 

= E A ^ J (\/^ - l )\J-\lA,v)] 2 +0 2 nch J (1) 

where cr 2 rich v = for unpunctured v (a^ represents the 

channel noise variance), <rL ch = for punctured v and 
{Ad„} is the edge perspective degree distribution of variable 
nodes. Similarly the EXIT function for the code component 
involving all check nodes is given by 

Ie,c(Ia,c) 

= 1 -Y,Pd,J(\/(dc - 1)[J- X (1 - ^,c)] 2 ) (2) 

d c 

where {pd c } is the edge perspective degree distribution of 
check nodes. 

IV. Semi-Structured E 2 RC-like Code Design 

In this section, we propose our design method for semi- 
structured E 2 RC-\\ke codes using EXIT charts. We consider 
the unstructured part and the structured part of E 2 RC codes 
shown in Fig. Q] as the two constituent code components. 
This code division for EXIT chart analysis is justified by the 
random edge interleaver between the two code components. 

We denote the set of variable nodes in the Tanner graph by 
V = Vi U V2 where V\ is the subset of nodes in H\ and V2 is 
the subset of nodes in H\2- The check node set is denoted by C. 
In the semi-structured E 2 RC codes, H2 has a base protograph 
structure of the form proposed in [6]. The base protograph 
shall be parameterized by the number of check nodes in 
it, denoted by M. The H2 part of semi-structured E 2 RC 
codes is obtained by simply replicating the base protograph 
an appropriate number of times. For example, the case of 
M = 8 is shown in Fig. Q] Check nodes are connected to 
the set V\ by a random interleaver (denoted II in Fig. [TJ. We 
shall frequently need to refer to the protograph representation 
of i?2- Let V p and C p denote the variable node set and the 
check node set in the protograph representation of f?2- Let 
e(vi) and e(c,) denote the set of edges connected to Vi € V p 
and Cj € C p respectively. We shall use £i,(cj) to denote the 
set of edges connecting Cj and the random edge interleaver 
and use £n(ci) to denote the set of edges connecting c, and 
V p , i.e., e(ci) = £z,(cj) Uen(ci). Given a protograph structure 
on H2, the problem of code design becomes one of finding 
good degree distributions for the variable nodes in V\ and 
that for the edges in yJ Ci eC v £L,{ci) (henceforth referred to 
as the left check degree distribution). In our examples, we 
only consider concentrated or near-concentrated total check 
degrees. We have found experimentally that these tend to give 
the best performance. Note that since the H2 part is fixed, this 
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implies that the left check degree distribution is also more or 
less fixed. Accordingly in our design process we experiment 
with a few check degree distributions and focus on optimizing 
the degree distribution for the nodes in V\. 

We explain our design method in the context of the binary- 
input AWGN (BIAWGN) channel. It can be adapted to the 
BEC and other channels in a straightforward manner. Suppose 
that we are given a channel noise variance <r^, the protograph 
specifying H2 and the left check degree distribution. The code 
design problem is to find the degree distribution {A<j„ , v £ V\} 
so as to minimize the gap between code rate R and channel 
capacity C (corresponding to cr%)> while constraining the 
maximum variable node degree to be c?„ max . Denote the EXIT 
function of the structured part by Ie,s(Ia,s)- For a given 
{\dv,v £ Vi}, the EXIT function of the unstructured part 
(see Fig. [TJ can be expressed as (according to (HJ) 

lE,uns(lA,unS, &mch) 

= E X ^ J (\/W' - W'^uns)] 2 + < ch ) 0) 
d v 

The code design or optimization problem is formulated as 
minimize : C — R 

d v , max 

subject to : 1. A<j„ = l,Ad„ > 

d„=l 

2- lE,uns{lA,uns) > Ia,s(Ie,s) 

for I A ,unS = Ie.s e [0, 1) 

Here, the second constraint is the zero-error constraint by 
ensuring the tunnel between the two EXIT curves. It is easy 
to see that minimizing C — R for a fixed a 2 n is equivalent 
to maximizing IT^ ^ or v ^ ^ ^he com P uta - 

tion of Ie,s(Ia,s) will be elaborated on in Section IIV-AI 
lE,uns(lA,uns) is a linear function of {X dv ,v £ as in ©. 
Therefore by a fine enough discretization of the interval [0, 1), 
we can express the above optimization as a linear program. 

In practice, to set up the second constraint, we need to 
find the inverse map Ia,s(Ie.s) by using linear interpolation. 
We have found that a large number (say 10 4 ) of {Ia,s, Ie.s) 
pairs for the function Ie,s(Ia,s) are necessary to ensure the 
accuracy of the inverse map Ia,s(Ie,s) an d the solution to 
the optimization problem. By solving the above optimization 
problem at a certain channel parameter a 2 , we get a code of 
rate corresponding to the {A^„} returned from the optimiza- 
tion. To get an optimized code at rate R , we need to solve 
the above optimization problem at closely spaced channel 
parameter levels below the Shannon limit corresponding to R Q 
until we get a code rate close enough to R Q . This necessitates 
numerous computations of Ie,s(Ia,s) an d motivates the need 
for a fast method for computing Ie,s{Ia,s)- 

A. A New Method for Computing EXIT Function of the 
Structured Part 

The usual approach for finding the EXIT function of a 
constituent code component is proposed in [20] by using 
Monte Carlo simulations. A large number of Monte Carlo 



simulations are needed for obtaining smooth EXIT functions. 
Moreover, this needs to be repeated at many different channel 
parameters. This makes the process rather time-consuming. 

Here, we present a fast and accurate method for computing 
EXIT functions of structured code components of LDPC 
codes, such as the structured part of E 2 RC codes and that 
of IRA codes, without resorting to Monte Carlo simulations. 

For convenience, we use the notation En = U t , ie v p e(i)i) 
and E L = U Cie c p £i(c l ). Note that \J Vi£Vp e{vi) = 
U Ci ec £r{ci)- Suppose that the a priori inputs carried on e £ 
El have average mutual information I am a nd that v £ V p has 
channel inputs parameterized by cr 2 nch v . We are interested in 
finding Ie, the average mutual information associated with the 
extrinsic outputs carried on e £ El after iterative decoding. 
For an edge e connected to node Wj(cj), we shall use the 
notation 7^* (likewise I^ e ) to denote the mutual information 
describing the a priori inputs on it and I E l (likewise 7g ) 
the mutual information describing the extrinsic outputs on it. 
We set up the following system of equations for the given 
structured code component. For e £ En and Vi £ V p , we have 




Similarly, for e £ En and Cj £ Cj 



I E\e = ^-J{{ E [J-^-Iam)? 

V V e'G£L(ci) 



e'ee R ( Cl )\{e} J J 

and for e € El and a £ C p , 




For each edge e £ En, there are two equations in the form 
of (O and (0 respectively and two unknown variables 1^ (or 
lA e ), Ie e ( or ^Ae) a ssociated with it; while for each edge 
e £ El, there is one equation in the form (O and one unknown 
variable I'g associated with it. We want to compute 1^ for 
e £ El- There are totally 2|£?r| + \El\ equations and the 
same number of unknown variables involved in this system of 
equations. Note that the specific expressions for this system 
of equations are totally dependent on the structure of the code 
component. The main idea behind our method for computing 
the EXIT function is to find the solution to this system of 
equations for a given value of I a. in a nd channel parameter. 
We now present an intuitive method for solving this system 
of equations, which works in an iterative manner by applying 
the sequence of updates described in equations (01, © and 
(O. The details are given below. 

1) Problem Instance. Given a structured code component, 
solve the system of equations described in ©, (O and (O 
above. The unknown variables involved in this system of 
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equations are T£ e (or I A 3 e ), -f^ e ( or f° r a ^ e e ^fl ar, d 
7g for all e 6 El. The known variables are 1% = Iaau for 
all e G El, and the channel parameter from which cr^t 
can be determined for each Uj. 

2) Initialization. Initialize all unknown variables to be 0. Set 
a small value of tthresh = 10 -6 . 

3) Iterative Updates. 

(a) Check node update. For all e G -Er, compute 7g using 
equation 10. Check to see whether the norm of the 
difference between this newly computed set of 7^ e and 
the previously computed ones is smaller than ethresh- If 
yes, then terminate; otherwise, set I A 3 e = /J? if and 
q are connected by e. 

(b) Variable node update. For all e £ -Erf, compute 7^f using 
equation ©. Set l^ 3 = F£ if Cj and are connected 
by e. Go to step 3(a). 

4) Compute I'g e for e G El. For all e G El, compute I'g e 
using equation (0. The average of these I'g e is denoted by 
Ie and (lA,in,lE) is a point on the EXIT function. 

The method can be adapted for computing EXIT functions 
over other channels by using appropriate update equations. 
Moreover, it can be used to compute the EXIT function of the 
structured part of other codes that have a succinct protograph 
representation such as IRA codes. 

We demonstrate the effectiveness of our method by compar- 
ing EXIT functions computed by our method and by the Monte 
Carlo simulation based method for two cases: the structured 
part of E 2 RC codes and that of IRA codes on BIAWGN 
channels respectively. The structured part of the E 2 RC code 
has a protograph structure of size 128. All check nodes have 
degree 8. To get smooth curves, we apply 10 6 a priori inputs 
in Monte Carlo simulations for computing each point on the 
curves. As shown in Table [I] the maximum absolute error 
(MAE) between the EXIT functions computed using the two 
methods is less than 0.0072. 

TABLE I 

Comparison of approaches for computing EXIT functions 



AWGN: noise variance— 0.95775;10 4 (Ia^e) pairs generated 




E' A B,C 


IRA 


Method 


Proposed 


Simulation 


Proposed 


Simulation 


Computing time (s) 


3.7 


24596 


0.6 


175886 


MAE 


0.0072 




0.00719 





B. Code Design Examples 

In our first example, we design a semi-structured E 2 RC 
code with d„ imax = 7 and check degree distribution of 
pe = 0.339623, p7 = 0.660377. Thus, from a complexity 
perspective these codes are comparable to the first design 
example in Section V in [10]. The mother code is of rate 0.5 
and the H2 part has a protograph structure of size M = 32. 
Our optimized code (referred to as code 0) is specified by A3 = 
0.4243, A7 = 0.5757 for v G V\ and has an asymptotic gap 
of 0.38 dB to capacity at rate 0.5. Fig. [3] gives the simulation 
results of this code of block length 2048 bits generated by 
the algorithms in [24] [25]. For comparison, we also list the 
simulation results of the reference code and original E 2 RC 
code from [10]. In this paper, our codes follow the designed 



puncturing patterns of original E 2 RC codes in [6] to get all 
puncturing code rates. From the simulation results it is clear 
that our code is better than the codes in [10] for all code 
rates. In particular they do not suffer from the high error floor 
problem of original E 2 RC codes at the mother code rate. 

In our second example, we design another semi-structured 
E 2 RC codes with concentrated check degree 8 and d v ^ nax = 
20. The optimized code (referred to as code 1) is given by 
A 3 = 0.305825, A 7 = 0.213474, A 8 = 0.181737, A 20 = 
0.298964 for v G Vi and it has an asymptotic gap of 0.217 
dB to capacity at rate 0.5 which is smaller than code 0. This 
is expected since d v ^ max is higher in this case. The simulation 
results of code 1 of block length 16384 bits are given in Fig. |4] 
Also given are the simulation results of the E 2 RC code that 
is constructed according to the degree distribution specified 
in [6] (referred to as original E 2 RC). It shows that our code 
achieves slightly better performance at rates near mother code 
rate but suffers a little at higher code rates. 

We use the following terminology in this paper. The pre- 
dicted threshold refers to the decoding threshold from asymp- 
totic code performance analysis and the measured threshold 
refers to the channel parameter where the code achieves BER 
= 10~ 4 in simulations. For our code of length 2048 bits, 
the measured threshold at rate 0.5 is 1.47 dB which is 0.9 dB 
away from the predicted one. For our code 1 of length 16384 
bits, this gap is only 0.4 dB. 



code 
- code in [9] 
original E 2 RC code 




3 3.5 4 4.5 
Eb/No (dB) 



\\\\ v 


•* code 
— 9 — code in [9] 
-A- original E 2 RC code 


> \ 


■ *v 

V - - - X. 



1.5 2 2.5 3 3.5 4 
Eb/No (dB) 



Fig. 3. Comparison between our code and the reference code in [10] of 
block length 2048 bits. The code rates are 0.5, 0.6, 0.7, 0.8 and 0.9 from left 
to right. The figure on the top (bottom) corresponds to BER (FER). 



C. Puncturing Performance Analysis and Joint Optimization 
of Semi-Structured E 2 RC Codes 

The puncturing performance of a given code is specified 
in terms of its decoding thresholds at all punctured code 
rates. The given semi-structured E 2 RC codes are specified 
by ^d v > v G V\ and the knowledge of the protograph structure 
of the structured part. For a given channel parameter, we can 
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compute the two EXIT functions using ^ and the approach 
of Section IIV-AI Note that when computing EXIT functions 
at different puncturing code rates, we follow the designed 
puncturing pattern of E 2 RC codes. The decoding threshold 
at a given code rate is determined by finding the channel 
parameter where the two EXIT curves (computed under the 
puncturing pattern at that rate) just begin to separate. 




0.5 0.7 0.9 1.1 1.3 1.5 1.7 1.9 2.1 2.3 2.5 2.7 2.9 3.1 3.3 3.5 3.7 3.9 



Eb/No (dB) 

Fig. 4. Comparison between our second code example and original 
E 2 RC code in [6] of block length 16384 bits. The code rates are 
0.5, 0.5714, 0.6667, 0.8 and 0.8889 from left to right. 

Our puncturing performance analysis of code 1 suggests that 
it has asymptotic decoding thresholds of around 0.40, 0.85, 
1.40, 2.45, 3.44 dB at rates ^, ^, j|, ^ and § respectively. 
The measured thresholds at these rates for the code of block 
length 16384 bits based on the simulation results in Fig. [4] 
are around 0.80, 1.22, 1.75, 2.78, 3.76 dB, which are consistent 
with the predicted ones with gaps uniformly around 0.35 dB. 
Joint Optimization of Semi-Structured E 2 RC Codes 
We now demonstrate that we can design our codes such that 
they have a small gap to capacity at all puncturing rates. 
This is because our puncturing pattern is deterministic and 
allows the determination of the asymptotic threshold at any 
puncturing rate for any given {A^, ,V G Vi}. Let R be a 
specified set of code rates where we want to optimize the code. 
Let a(g, Ri) denote the channel noise parameter that is at a gap 
of g from the channel parameter corresponding to the Shannon 
limit at rate Ri. Let Ia,s(Ie,s, o~(g, Ri)) denote the plot of 
Ia,s vs. Ie,s under the puncturing pattern corresponding 
to rate Ri, at the channel parameter a(g, Ri). The notation 
lA,uns(lA,unSi &(g, Ri)) will be used analogously. We can 
formulate the joint optimization problem as minimizing the 
maximum gap to capacity at all rates in R as follows. 
Joint optimization algorithm 

for g = grain '• graax 

Solve the following linear program optimization problem 

d VjTnax 

EAd„ 

rt 

.1 i Uv 

d v —l 

d v , max 

subject to:l. Xd v = 1, Xd v > 0, 

d v = l 

2. lE,uns{lA,unS,e{g,Ri)) > Ia,s(Ie,S, &(g, Ri)) 

for I A ,unS = Ie.s £ [0, 1), for all Ji,eR 
if mother code rate corresponding to {A^, } is acceptable 



break; return {A^} and g. 
endif 
endfor 

Note that the second set of constraints is the zero-error 
constraint by ensuring an iterative decoding tunnel for all 
EXIT charts at all the required code rates. We can obtain all the 
required EXIT functions relatively quickly using our approach 
outlined previously under the puncturing patterns for each Ri. 
To obtain optimized {Xd v }, we basically keep increasing g 
until we get the desired code rate. The code specified by A<j„ 
is guaranteed to have asymptotic performance gap to capacity 
no larger than g at all code rates in R. 

We designed a semi-structured E 2 RC code that was jointly 
optimized across the rate range ^ ~ I, where M = 32, 
all check nodes have degree 8 and rft,. ma x = 20. The code 
is specified by A 3 = 0.309090, A 6 = 0.278794, A 20 = 
0.412116. Fig. gives the simulation results for the code 
of block length 16384 bits (listed as code 2). Also plotted 
are the simulation results in Fig. |4] for code 1 and original 
E 2 RC code from Section HV-BI The puncturing performance 
analysis suggests that code 1 has asymptotic performance 
gaps around 0.21,0.32,0.34,0.41,0.405 dB to capacity at 
rates j^, j^, an d § respectively while code 2 has much 
more uniform gaps of around 0.29,0.30,0.25,0.29,0.295 dB 
to capacity at these rates. The simulation results in Fig. [5] 
also suggest uniformly better performance of code 2 compared 
to code 1 across the range of rates. Moreover, the measured 
thresholds of the two codes at all rates in simulations are in 
good agreement with the predicted thresholds. At all code 
rates, the gap between the measured threshold and the pre- 
dicted threshold is around 0.35 dB. Finally, we note that code 
2 achieves better or at least the same performance as original 
E 2 RC code at all rates. 

The methods described in this section apply more gen- 
erally to codes that have a structured component with a 
protograph representation, such as IRA codes. We applied 
our method to the design of IRA codes as well and ob- 
tained jointly optimized codes with performance gaps around 
0.27,0.30,0.195,0.27,0.29 dB to capacity at rates ^ 
j^, and | respectively. Though not presented here, the 
simulation results of the IRA code are almost identical to the 
jointly optimized E 2 RC code (code 2) discussed above. 

To the best of our knowledge, a joint optimization algorithm 
that minimizes the gap to capacity simultaneously across all 
code rates has not been considered previously in the literature. 
In [2], the authors found the optimal puncturing patterns for 
two optimized mother codes (referred to as Ha code 1 and Ha 
code 2) and also gave their asymptotic puncturing performance 
using Gaussian approximation based density evolution. In 
Fig. [6] we compare the asymptotic thresholds of their codes 
with our code families. Note that our codes have the same 
maximum variable node degree and slightly smaller average 
variable node degree compared to the codes in [2]. We observe 
that the gaps to capacity for our codes remain more or less 
the same across the range of rates, whereas the codes in [2] 
exhibit a larger gap to capacity at higher code rates. 
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Fig. 5. Comparison between two semi-structured E 2 RC codes optimized 
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(code 2) and original E 2 RO code of block length 16384 bits. The code rates 
are 0.5, 0.5714, 0.6667, 0.8 and 0.8889 from left to right. 
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V. Protograph E 2 RC codes construction 

In this section, we introduce the construction of a class of 
structured E 2 RC-\\ke codes based on protographs. The basic 
idea in these codes is to impose a protograph structure on 
the systematic part Hi of the parity-check matrix as well (in 
addition to the protograph structure on H2) . We obtain a family 
of protographs with asymptotic gaps to capacity no larger than 
0.28 dB across a wide range of rates when the maximum 
variable node degree is twenty. These codes have excellent 
finite length performance as well. In this part of the work, we 
use the reciprocal channel approximation of density evolution 
for computing the threshold for a given protograph [9] [19] [13]. 
The construction algorithm goes as follows. 

1) Find a good high-rate protograph. Using density evolution 
we first identify a high-rate protograph (starting protograph) 
with a low threshold. 

2) Use check-splitting to obtain lower-rate protographs with 
low thresholds. From the starting protograph, we perform 
check splitting in a systematic manner to obtain a family of 
good rate-compatible punctured protographs where the parity 
part has the E 2 RC structure. These correspond to different 
code rates of the rate-compatible code family. 

3) Construct the LDPC code by replacing protograph edges 
with carefully chosen circulant permutations. With the pro- 
tograph of mother code rate constructed from above steps, a 
larger graph defining the LDPC code is constructed by replac- 
ing the protograph edges with appropriately chosen circulant 
permutations by using techniques in [26] [24] [25] [27]. 



In the sequel we shall attempt to explain the construction 
process by means of an example. However, it should be clear 
that the techniques are applicable in general. 

A. Starting Protograph 

Let the desirable code rate range for the code family be 
Rmin < R < Rmax ■ A high-rate protograph of size Mq x Nq 
(Mo - number of check nodes, No - number of variable 
nodes) with low threshold serves as the starting protograph. 
The mother code protograph is of size M x N such that 
^ir 1 < R ™n, > Rmax and N - M = N Q - M . 

These conditions guarantee that the desirable code rate range 
is achievable by the construction. In addition, these parameters 
should be kept relatively small (less than 50) to keep the 
construction complexity manageable. We impose the constraint 
that the degree of the variable nodes in the starting protograph 
is at least three to avoid a high error floor. We perform an 
exhaustive search using the reciprocal channel approximation 
of density evolution [19] [13] to find a protograph with low 
threshold of size Mo x iVo and with variable node degree 
between 3 and a maximum degree d v ^ max . In our example, 
Rmin = 0.5, Rmax = §, d v , max = 20 and the size parameters 
are decided as M = 1,N = 9 and M = 8,N = 16. Our 
example starting protograph consists of one check node and 
nine variable nodes of degree {20,8,3,3,3,3,3,3,3} respectively, 
i.e. each variable node is connected to the single check node 
by multiple edges. It has a threshold (computed according to 
[19][13]) of 3.27 dB which is 0.24 dB away from capacity. 

B. Check-Splitting 

The operation of check-splitting (also used in [28]) on a 
check node c of degree d in the protograph G\ proceeds as 
shown in Fig. [7] We split c into two new check nodes c\, of 
degree d\, and C2, of degree d2, such that d = d%+d2- Next we 
introduce a new variable node v§ and introduce edges ci — Ve 
and C2 — Ve, so that the degree of vq is two. The resultant 
protograph G3 is of lower rate than G\. 

G, G, 




Fig. 7. Check node c with degree-5 is divided into ci and C2. The new 
check nodes c\ and C2 are connected by a new degree-2 variable node Ug. 



Starting with the high-rate starting protograph, we apply 
check-splitting repeatedly in a specific manner, and finally 
arrive at the protograph of low rate mother code. In fact, the 
protographs produced in the check-splitting process form a 
family of rate-compatible protographs if we consider the newly 
added degree-two variable nodes in check-splitting as parity 
nodes providing incremental redundancy. However, the check- 
splitting needs to be done carefully in order to have good code 
performance across all rates. We note that performing density 
evolution on G\ in Fig. [7] can predict the threshold of G3 
under puncturing (upon puncturing Vq, the resultant G2 has 
the same asymptotic decoding threshold as G\) [29]. 
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C. Constructing Protographs with E 2 RC-like Structure 

We shall call the original variable nodes of degree at least 
three in the starting protograph, old nodes and the variable 
nodes of degree two introduced in check-splitting, new nodes. 
For a given check node, we define its old (new) node degree to 
be the number of connections to the old (new) nodes. When we 
split Co into coi and C02, a decision needs to be made on how 
the connections of cq are divided between them. To obtain the 
E 2 RC structure, coi is allocated all of co's new node degree; 
the parity node newly introduced in check-splitting has one 
connection to both coi and cq2 (for a proof see the Appendix). 
The old node degree also needs to be divided between cqi 
and C02 in a manner that ensures that the threshold of the new 
protograph is low. We discuss it in more detail in Section [V-DI 

At each stage of the construction, we perform check- 
splitting on all check nodes in the current protograph. We 
use the example starting protograph from Section IV-AI to 
demonstrate the process. Here Mq = 1 and M = 8, so we 
shall have log 2 8 = 3 construction stages. Let M„ s denote the 
number of check nodes in the protograph at the beginning of 
stage n s . In stage n s , we perform check-splitting on all M ns 
check nodes so that a set of M ns protographs of decreasing 
rates are generated at this stage. The order in which the check 
nodes are split can affect the thresholds at those rates. 

We now show the first and second splitting stages for our 
starting protograph. In the first stage, check-splitting on the 
single check node generates a new protograph of rate y^, 
shown in Table [TTJ In the second stage, there are two check 
nodes in the protograph. Density evolution analysis tells us 
that performing check-splitting on cqi first gives a protograph 
of rate yy with a better decoding threshold. So in this stage, 
we first split cqi to generate a protograph of rate yy and then 
split C02 to generate a protograph of rate y| (see Table ITTTb - 
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The third stage proceeds in a similar manner. Refer to Fig. 
[8] for a graphical illustration of the construction process. Note 
that at the end of stage n s of the algorithm, the M„ s newly 
added parity nodes are all 1-SR. Also a k-SR node at the end 
of stage n s — 1 would become (k+l)-SR at the end of stage 
n s (as shown in the appendix). In this way the construction 
ensures that half of the parity nodes in the final protograph 
are 1-SR, one-fourth are 2-SR and so on, i.e., our construction 
results in protographs with E 2 RC structure. The parity nodes 
of the resultant mother code protograph are punctured in the 



inverse order in which they were added in the construction in 
order to obtain higher puncturing rates. 




Fig. 8. Protograph construction from check-splitting. 



D. Deciding Splitting Patterns 

In check-splitting, the connections between a given check 
node and the old nodes can be divided among the two new 
check nodes in many ways. Let Cq be the check node to be 
split. Let so be the vector of connections between Co and 
the old nodes, e.g., so = [20 8 3 3 3 3 3 3 3] in our 
example. The splitting pattern refers to the set of vectors 
soi = [10 4 2 1 2 1 2 1 2] and s 02 = [10 4 1 2 1 2 1 2 1] 
(shown in Table [TTJ that determine the connections between 
Cox an d c 02 and the old nodes. Thus so = «oi + s 02- It is 
easily seen that the number of possible splitting patterns are 
huge in each check-split and it's impossible to evaluate all of 
them in our construction. For example, when splitting sq into 
s i and s 02 , we have totally 20 1 x 8 1 x 3 7 = 349920 possible 
splitting patterns. Our search for good splitting patterns for 
each check-splitting is guided by two main points. 

a) Trade-off between the performance of high-rate and low- 
rate protographs. We have found that there exists a tradeoff 
between the performance of high-rate and low-rate protographs 
during the construction. For example, very low thresholds for 
the higher rate protographs typically come at the expense of 
higher thresholds for the low-rate protographs. 

b) Equal splitting patterns give good performance. Note that 
considering all possible splitting patterns at all possible stages 
is essentially computationally infeasible. We have found that 
splitting patterns that split the connections roughly equally 
between the two new check nodes in each check-splitting 
result in good thresholds across all code rates in the family. 
This reduces the search space a lot and it becomes possible 
to perform density evolution analysis to determine proper 
splitting patterns at each stage. 

In Table IIV1 we present one of the mother code protographs 
that we have constructed from the example starting protograph 
from Section IV-AI using the construction algorithm described 
above. Also shown is the gap to Shannon limit for the 
protograph at different puncturing rates | ~ y^. The gap to 
capacity remains between 0.235-0.278 dB across the range 
of rates. In fact, we could choose other splitting patterns such 
that the protograph family has an even lower threshold at some 
high code rate during the early stage of the construction, but 
it would come at the expense of higher thresholds for the 
lower code rates in the later stage of the construction. These 
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protographs are available at [30] (Due to lack of space we are 
unable to include them here). 

TABLE IV 
Mother code Protograph-1 
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E. Results 

We compared the asymptotic thresholds of the code family 
represented by protograph-1 above and the AR4JA code family 
in [28]. For the three common code rates i, |, and |, our 
protograph is better than the AR4JA family with performance 
gains of 0.17,0.143 and 0.12 dB respectively. The average 
variable node degree of our codes is a little higher than that 
of AR4JA family. Note however that the codes in [28] are not 
rate-compatible punctured codes. 

We constructed protograph E 2 RC codes with block length 
16384 bits by replacing each edge in protograph-1 by an 
appropriate circulant permutation using the algorithms pro- 
posed in [24] [25]. Fig. [9] gives the simulation results of this 
protograph E 2 RC code, the jointly optimized semi-structured 
E 2 RC code (code 2) and original E 2 RC code of the same 
block length from Section IIV-CI From the asymptotic perfor- 
mance analysis for protograph E 2 RC code as shown on the 
bottom of Table IPVl and for code 2 in Section |TV-CI we see that 
protograph E 2 RC code is quite competitive to the optimized 
semi-structured E 2 RC code. The simulation results shown 
in Fig. [9] are consistent with the analysis. Moreover, both 
protograph E 2 RC code and optimized semi-structured E 2 RC 
code achieve better performance than the original E 2 RC code. 




0.5 0.7 0.9 1.1 1.3 1.5 1.7 1.9 2.1 2.3 2.5 2.7 2.9 3.1 3.3 3.5 3.7 3.9 



Eb/No (dB) 

Fig. 9. Comparison between protograph E 2 RC code, optimized semi- 
structured E 2 RC code and the original E 2 RC code of block length 16384 
bits at rates 0.5, 0.5714, 0.6667, 0.8 and 0.8889 from left to right. 



VI. Conclusions 

The E 2 RC codes were proposed in [6] as a promising class 
of rate-compatible codes. In this work we introduced semi- 
structured i? 2 i?C-like codes and protograph E 2 RC codes. We 
developed EXIT chart based methods for the design of semi- 
structured i? 2 i?C-like codes that allow us to determine near- 
optimal degree distributions for the systematic part of the code 
while taking into account the structure of the deterministic 
parity part. We presented a novel method for finding EXIT 
functions for structured code components that have a succinct 
protograph representation that is applicable in other scenarios 
as well. This allows us to analyze the puncturing performance 
of these codes and obtain codes that are better than the 
original construction. Using our approach we are able to 
jointly optimize the code performance across the range of rates 
for our rate-compatible punctured codes. Finally we consider 
E 2 RC-like codes that have a protograph structure (called 
protograph E 2 RC codes) and propose design rules for rate- 
compatible protographs with low thresholds. These codes are 
useful in applications since the protograph structure facilitates 
implementation. For both the semi-structured and protograph 
E 2 RC families we obtain codes with small gaps to capacity 
across the range of rates. 
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Appendix 

We first recall the precise construction rule. Let Co be a 
check node with a certain number of connections to new nodes 
of degree-2. When co is split into coi and C02, coi inherits all 
of co's connections to new nodes and both coi and C02 have 
one connection with the newly introduced variable node. 

Suppose that the starting protograph is of size Mo x Nq. 
At the end of stage k of the construction algorithm, it is clear 
that there will be M ■ 2 k check nodes and M (2 k - 1) new 
degree-2 parity nodes. Our aim is to show that this construction 
algorithm results in an H 2 part with the E 2 RC structure, i.e. 
at the end of stage k of the algorithm, there are Mq ■ 2 k ~ 1 
1-SR nodes, M ■ 2^ 2 2-SR nodes, . . . , and M k-SR nodes. 
We proceed by induction. 

Base Case. At the end of the first stage, we will have 2Mq 
check nodes and Mq new degree-2, 1-SR parity nodes. 

Inductive Step. Suppose that the statement is true at the end 
of stage k. We will show that it is true at the end of stage k + 1. 
To see this note that at the end of stage k + 1 we will have 
Mq ■ 2 k+1 check nodes formed by splitting the check nodes at 
the end of stage k. According to the construction algorithm, 
MQ-2 k check nodes will inherit the previous connections while 
the remaining will have just one connection to the Mq ■ 2 k 
newly introduced degree-2 variable nodes. Therefore we will 
have at least Mq ■ 2 k 1-SR nodes at the end of stage k + 1. 
Next we note that any node that was of type a-SR at the 
end of stage k will now become of type (q+1)-SR. This is 
because each of the check nodes it is connected to will have 
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one additional connection. This implies that at the end of stage 
k + 1 there will be M • 2 k 1-SR nodes, M • 2 fe ~ 1 2-SR nodes, 
. . . , and M (k+l)-SR nodes. This shows the required result. 

References 

[1] R. G. Gallager, Low Density Parity Check Codes. MIT press, 
Cambridge, MA, 1963. 

[2] J. Ha, J. Kim, and S. W. McLaughlin, "Rate-Compatible Puncturing of 
Low-Density Parity-Check Codes," IEEE Trans. Inform Theory, vol. 50, 
no. 11, pp. 2824 - 2836, Nov. 2004. 

[3] , "Puncturing for Finite Length Low-Density Parity-Check Codes," 

IEEE Int. Symp. Inform. Theory, Chicago, p. 152, Jun. 2004. 

[4] J. Ha. J. Kim, D. Kline, and S. W. McLaughlin, "Rate-Compatible 
Punctured Low-Density Parity-Check Codes with Short Block Lengths," 
IEEE Trans. Inform. Theory, vol. 52, no. 2, pp. 728 - 738, Feb. 2006. 

[5] G. Yue, X. Wang, and M. Madihian, "Design of Rate-Compatible 
Irregular Repeat Accumulate Codes," IEEE Trans. Commun., vol. 55, 
no. 6, pp. 1153 - 1163, Jun. 2007. 

[6] J. Kim, A. Ramamoorthy, and S. W. McLaughlin, "Design of Efficiently- 
Encodable Rate-Compatible Irregular LDPC Codes," IEEE Trans, on 
Communications, vol. 57, no. 2, pp. 365 - 375, Feb. 2009. 

[7] M. R. Yazdani and A. H. Banihashemi, "On Construction of Rate- 
Compatible Low-Density Parity-Check Codes," IEEE Comm. Letters, 
vol. 8, no. 3, pp. 159 - 161, Mar. 2004. 

[8] J. Kim, W. Hur, A. Ramamoorthy, and S. McLaughlin, "Design of Rate- 
Compatible Irregular LDPC Codes for Incremental Redundancy Hybrid 
ARQ Systems," IEEE Int. Symp. Inform. Theory, pp. 1139 - 1143, Jul. 
2006. 

[9] J. Thorpe, "Low Density Parity Check (LDPC) Codes Constructed from 
Protographs," JPL INP Progress Report 42-154, Aug. 2003. 

[10] S. Song, D. Hwang, S. Seo, and J. Ha, "Linear-Time Encodable Rate- 
Compatible Punctured LDPC Codes with Low Error Floors," IEEE VTC, 
pp. 749 - 753, May. 2008. 

[11] H. Jin, A. Khandekar, and R. McEliece, "Irregular Repeat-Accumulate 
Codes," Proc. 2nd Int. Symp. Turbo Codes and Related Topics, Brest, 
France, pp. 1 - 8, Sep. 2000. 

[12] A. Roumy, S. Duemghar, G. Caire, and S. Verdu, "Design Methods for 
Irregular Repeat-Accumulate Codes," IEEE Trans. Inform. Theory, vol. 

50, no. 8, pp. 1711 - 1727, Aug. 2004. 

[13] T. Richardson, "Multi-Edge Type LDPC Codes," presented at the 
Workshop honoring Prof. Bob McEliece on his 60th birthday, California 
Institute of Technology, Pasadena, California, May. 2002. 

[14] T. J. Richardson, M. A. Shokrollahi, and R. L. Urbanke, "Design 
of Capacity-Approaching Irregular Low-Density Parity-Check Codes," 
IEEE Trans. Inform. Theory, vol. 47, no. 2, pp. 619 — 637, Feb. 2001. 

[15] T. J. Richardson and R. L. Urbanke, "The Capacity of Low-Density 
Parity Check Codes Under Message-Passing Decoding," IEEE Trans. 
Inform. Theory, vol. 47, no. 2, pp. 599 - 618, Feb. 2001. 

[16] M. G. Luby, M. Mitzenmacher, M. A. Shokrollahi, D. A. Spielman, and 
V. Stemann, "Practical Loss-Resilient Codes," Proc. 29th Annu. ACM 
Symp. Theory of Computing, pp. 150 - 159, 1997. 

[17] M. A. Shokrollahi, "New Sequences of Linear Time Erasure Codes 
approaching the Channel Capacity," Proc. 13th International Symposium 
Applied Algebra, Algebraic Algorithms, and Error-Correcting Codes, pp. 
65 - 76, 1999. 

[18] S. ten Brink and G. Kramer, "Design of Repeat-Accumulate Codes for 
Iterative Detection and Decoding," IEEE Trans. Signal Processing, vol. 

51, no. 11, pp. 2764 - 2772, Nov. 2003. 

[19] S.-Y. Chung, On the Construction of Some Capacity-Approaching Cod- 
ing Schemes. Ph.D. dissertation, Massachusetts Institute of Technology, 
Cambridge, Massachusetts, Sep. 2000. 

[20] S. ten Brink, "Convergence Behavior of Iteratively Decoded Parallel 
Concatenated Codes," IEEE Trans. Commun., vol. 49, no. 10, pp. 1727 
- 1737, Oct. 2001. 

[21] S. ten Brink, G. Kramer, and A. Ashikhmin, "Design of Low-Density 
Parity-Check Codes for Modulation and Detection," IEEE Trans. Com- 
mun., vol. 52, no. 4, pp. 670 - 678, Apr. 2004. 

[22] A. Ashikhmin, G. Kramer, and S. ten Brink, "Extrinsic Information 
Transfer Functions: Model and Erasure Channel Properties," IEEE 
Trans. Inform. Theory, vol. 50, no. 11, pp. 2657 - 2673, Nov. 2004. 

[23] E. Sharon, A. Ashikhmin, and S. Litsyn, "Analysis of Low-Density 
Parity-Check Codes Based on EXIT Functions," IEEE Trans. Commu- 
nications, vol. 54, no. 8, pp. 1407 - 1414, Aug. 2006. 



[24] T. Tian, C. R. Jones, J. D. Villasenor, and R. D. Wesel, "Selective 
Avoidance of Cycles in Irregular LDPC Code Construction," IEEE 
Trans. Commun., vol. 52, no. 8, pp. 1242 - 1247, Aug. 2004. 

[25] A. Ramamoorthy and R. Wesel, "Construction of Short Block Length 
Irregular Low-Density Parity-Check Codes," IEEE International Con- 
ference on Communications, vol. 1, pp. 410 - 414, 2004. 

[26] X.-Y. Hu, E. Eleftheriou, and D.-M. Arnold, "Progressive Edge-Growth 
Tanner Graphs," IEEE GlobeCom, vol. 2, pp. 995 - 1001, Nov. 2001. 

[27] W.-Y. Weng, A. Ramamoorthy, and R. Wesel, "Lowering the error 
floors of irregular high-rate LDPC codes by graph conditioning," IEEE 
Vehicular Technology Conference, vol. 4, pp. 2549 - 2553, Sep. 2004. 

[28] D. Divsalar, S. Dolinar, and C. Jones, "Construction of Protograph LDPC 
Codes with Linear Minimum Distance," IEEE Int. Symp. Inform. Theory, 
pp. 664 - 668, Jul. 2006. 

[29] H. Pishro-Nik and F. Fekri, "Results on Punctured Low-Density Parity- 
Check Codes and Improved Iterative Decoding Techniques," IEEE 
Trans. Inform. Theory, vol. 53, no. 2, pp. 599 - 614, Feb. 2007. 

[30] http://www.ece.iastate.edu/~cshi/web.pdf. 



10 



